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Abstract 

We investigate the impact parameter dependent parton distributions for a relativistic composite 
system in light-front framework. We take an effective two-body spin-1/2 state, namely an electron 
dressed with a photon in QED. We express the impact parameter dependent parton distributions in 
terms of overlaps of light-cone wave functions. We obtain the scale dependence of both fermion and 
gauge boson distributions and show the distortion of the pdfs in the transverse space for transverse 
polarization of the state at one loop level. 
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I. INTRODUCTION 

Impact parameter dependent parton distributions q(x, b L ) have been introduced re- 
cently as a physical interpretation of generalized parton distributions (GPDs) [3] in terms of 
probability densities in the impact parameter space. It is known that GPDs are off-forward 
matrix elements of light-cone bilocal operators and they do no have a probabilistic interpre- 
tation like ordinary parton distributions. q(x, b- 1 ) can be expressed as Fourier transform of 
GPDs with respect to the transverse momentum transfer (when the longitudinal momentum 
transfer is zero) and they give the distribution of partons in the transverse position space. 
In fact they obey certain positivity constraints and thus it is legitimate to associate the 
physical interpretation as a probability density. This interpretation is not limited by rela- 
tivistic effects in the infinite momentum frame. Another interesting aspect is that when the 
state is polarized in the transverse direction, the unpolarized impact parameter dependent 
pdf is distorted in the transverse plane This distortion can be linked to the transverse 
single spin asymmetries observed in transversely polarized Lambda production in pp or pp 
collisions. Recently it has been shown in the framework of the scalar diquark model that the 
impact parameter space asymmetry together with the final state interaction of the active 
quark gives rise to the Sivers effect |4| in momentum space Q]. 

Impact parameter dependent pdfs have been investigated in chiral quark model for the 
pion p, in the constituent quark model |7fl, in terms of a power-law ansatz of the light- 
cone wave function for the pion sj], in the context of investigating the color transpar ency 



phenomena |9[, in the transverse lattice framework for the pion |10| and on the lattice 11]. 
In this work, we investigate both fermion and gauge boson distributions for a relativistic 
composite system in the light-front framework, taking into account the correlation between 
different Fock components of the light-cone (or light-front) wave function in light-front gauge. 
It is well known that the light-front formalism is especially suitable to investigate many 
aspects of relativistic bound states because of the Galilean transverse boost operators and 
the triviality of the vacuum when a cutoff is imposed on the longitudinal momentum. 
We take an effective spin | system, namely an electron dressed with a photon in QED. 
It is known that such a model is self consistent and has been used to investigate the helicity 
structure of a composite relativistic system Q|. The state can be expanded in Fock space 
in terms of light-cone wave functions. The light-cone wave functions in this case can be 
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obtained from perturbation theory, and thus their correlations are known at a certain order 
in the coupling constant. However, their general form provides a template for parametrizing 
the more realistic composite system, namely the hadronic wave function Instead of an 
electron, one can also consider a state like a dressed quark in perturbative QCD. In fact this 
can be thought of as a field theoretic parton model where the partons, or the quarks and 
gluons, have mass, intrinsic transverse momenta and they interact. Previous studies have 
shown that this gives an intuitive picture of the DIS structure functions and scaling violations 
[15I ] and is suitable to address issues related to the spin and orbital angular momentum of 



the nucleon 







171 ]. Such a state has also been used to investigate the twist three GPDs 



in terms of overlaps of light-cone wave functions 18(. The aim of the present work is to 
extend these studies to the impact parameter dependent fermion and gauge boson pdfs in 
order to obtain a qualitative behavior in the impact parameter space within the framework 
of perturbation theory. 

The plan of the paper is as follows. In section II, we give the definition of the impact 
parameter dependent pdf, In sections III and IV we calculate the fermion and gauge boson 
distributions respectively. Summary and conclusions are given in section V. 



II. IMPACT PARAMETER DEPENDENT PARTON DISTRIBUTIONS 

Impact parameter dependent pdfs are defined by considering a transversely localized state 

Q: 

I P+, & = 0, A) = N J |^ I P + , P 1 , A), (2.1) 

where | P + , P -1 , A) is the light cone helicity eigenstate with momentum P and helicity A, 
N is a normalization factor. For a state with total momentum P + the transverse center of 
momentum R 1 - is defined as, 

R 1 = —— f c/x-rf 2 x ± 6 ++ x ± , (2.2) 



2P+ 

where is the energy-momentum tensor. 

The states | P + , P^, A) can be thought of as helicity eigenstates in the infinite momentum 
frame [jj. Due to the fact that the light-front transverse boost B 1 - behaves similar to the 
non-relativistic Galilean Boost generators, it can be shown that the states | P + , R 1 - = 0, A) 



defined above are simultaneous eigenstates of the longitudinal light cone momentum P + , 
transverse position operator R 1 - and light cone helicity J z . The transverse position operator 
P- 1 in fact is related to B 1 : P- 1 = -^+B L . 

The impact parameter dependent pdf is defined as, 

q(x, b L ) = (P + , R 1 = X , A | O q (x, 6 X ) | P + , R ± = X , A) (2.3) 

with 

O q (xX) = J ^{- V Y M)l + ^{\M)e^ P+y ^ (2.4) 

b ± is the impact parameter, which is the transverse distance of the active quark from the 
center of mass. We have taken the light-front gauge, A + = 0. Instead of the fermion 
operator, one can also have a gauge boson operator 

O g (x, b±) = / %F +V {~\, b^)F\( y T , b^ p+ y~ (2.5) 

to define the gauge boson pdf g(x, b^) similar to Eq. (|2.3|) . 

It can be shown that q(x, b- 1 ) can be expressed as a Fourier transform of the GPD 
H g (x,0,A 2 ) [l| : 

q(x, b±) = H q (x, b L ) = J ^ e - lb± - A± H g (x, 0, A 2 ). (2.6) 

where A 2 is the total momentum transfer and it is purely transverse. One can obtain a 
similar relation for the gauge boson counterpart. 



III. FERMION DISTRIBUTIONS 

We first calculate the off-forward matrix elements, parametrized in terms of the twist-two 
GPSs H q (x,A 2 ) and E q (x,A 2 ) in the standard way [21J and using the light cone spinors 
0- We have, 

/ tr eH "" (p ' Tl ^-Y h+ ^Y ] 1 p T) = H *( x > A2)) (3 - 1} 

where H(x, A 2 ) = H(x, 0, A 2 ). The momentum of the initial state is P M and that of the final 
state is P lfl . The average momentum between initial and final state is then P M = EH±EJL_ 
The momentum transfer is given by A M = P'^ — P M , P' ± = —P± = 4p. We take the 



skewedness £ = 0, in other words, the momentum transfer is purely transverse. We exclude 
the impact parameter space representation of GPDs for non-zero skewedness ^9| from our 
analysis here. 

We take the state | P, a) to be a dressed electron consisting of bare states of an electron 
and an electron plus a photon : 



P,a) =JV[6 t (P,<r) |0> 



dk+dtki / 2(27r) 3p +5 3 ( p _ h _ h) 



(j> 2 (P, a | k\, <ti; k 2 , A 2 )6 f (fci, c7i)a t (fc 2 , A 2 ) | 0) 



(3.2) 



Here and are bare photon and electron creation operators respectively and 2 is the 
multiparton wave function. It is the probability amplitude to find one electron plus photon 
inside the dressed electron state. 

We introduce Jacobi momenta x iy such that Y^i^i — 1 and J2i ( li ± = 0. They are 
defined as 

Also, we introduce the wave function, 



Qi 



kf - x i P ± . 



(3.3) 



(3.4) 



which is independent of the total transverse momentum P 1 - of the state and boost invariant. 
The state is normalized as, 

(P', A' | P, A) = 2(2tt) 3 P + (5 Ai v(5(P + - P' + )5 2 (P ± - P !± ). (3.5) 

The two particle wave function depends on the helicities of the electron and photon. Using 
the eigenvalue equation for the light-cone Hamiltonian, this can be written as Q], 



x(l — x) 



1 



(^) 2 + - x) 2 ^2(2^)3 
. _ i (1 — x 



X(T\ 



Q 



x 



x 



x 



(3.6) 



m is the bare mass of the electron, a 2 = — er 1 and er 1 = a 2 . In our case, A 2 = — (A -1 ) 2 as 
£ = 0. J\f gives the normalization of the state 



N 2 



a 
2^ 



1-e 



dx 



l + x 2 Q 2 l+x 2 

l °9— ~ ~ ~ + I 1 -x) 



x 



X 



(3.7) 



within order a. Here e is a small cutoff on x. We have taken the cutoff of the transverse 
momenta to be Q 2 , a 2 is a scale which we have taken to be m 2 << a 2 << A 2 , it mimics the 
factorization scale in QCD separating hard and soft dynamics |15|]. The above expression is 
derived using Eqs (|3.5|) . (|3.2j) and (|3.6|) . 

Contribution to H q (x, A 2 ) comes from one particle and two-particle sectors. The one- 
particle sector contributes only when x = 1 and is a delta function. This receives correction 
upto order a from the normalization condition of the state. The contribution from the 



two-particle sector can be written as an overlap of wave functions 

A/" 2 E / d'q^ix, + (1 - x)A ± )^, ^ 
Using the expression Eq. ()3.6|) we get 



18, 



22] 



(3.8) 



H q (x,A 2 ) 



+ 



N 2 5(l - x) 



X ^(q L f + m 2 (l-xf 



(2tt) 3 
+ (l + a: 2 ) J d 2 q J 

One can write this as, 

H q (x,A 2 ) = Af 2 {5(l- 



d 2 q~ 



((q^) 2 + m 2 (l - x) 2 )((g^ + (1 - x)A ± ) 2 + m 2 (l - x) 2 ) 



q^.A J 



x 



—m 2 (l — x) 2 



d 2 q J 



((q^ 2 + m 2 (l - x) 2 )((g x + (1 - x)A ± ) 2 + m 2 (l - x) 2 ) 

d 2 q L 
L 2 

fq 1 - 1 + x 2 r r d 2 q ± 
X) J TJT 2 + 2(1- x)U 
d 2 q~ 



.(3.9) 



(2vr) 3 
d 2 q ± - 



L\L 2 
l-x) 2 (A L ) 2 



+ m 2 (l - a;) 3 



Li 



where L\ = (q L ) 2 + m 2 (l 



x 



and L 2 = (q A 



L\L 2 - _ 
1 -x)A ± ) 2 + m 2 (l 



(3.10) 
x) 2 . In the forward 



limit using the normalization condition of the state we get, 



H q (x,0) = 5(1 -x) + 



-5(1 -x) 



a 

27T* 



1 



x 2 1 Q 2 

- x jl z 



1 + x 1 



+ 1 - X 



X 



d <T^ l ° 9 ^ 



i + y 2 
1-3/ 



(3.11) 



Integrating over x, one gets Jq 1 dxH q (x, 0) = -Fi(O) = 1, where Fi(0) is the form factor at 
zero momentum transfer. Replacing a by a s Cf and neglecting the mass, eq. (|3.11J) reduces 
to the quark distribution function of a dressed quark and the coefficient of the logarithmic 
term gives the splitting function P qq At non-zero A^, in the limit x — > 1, we have 

-V O 2 r 1 -J- 3 

(3.12) 



# 9 (£,A 2 ) 



1+^ 2 , 3 , 



The plus '+' prescription is defined in the usual way. There is no singularity at x — ► 1 
and H q (x,A 2 ) is independent of A 1 - at x — > 1. The impact parameter dependent parton 
distribution is obtained by taking a Fourier transform with respect to A -1 . At x — > 1 , we 

get 

n q (xX) = j '^e~ ib± -* X H g (x,A 2 ) 

= 6^){6(1 -x) + ^log^ [1±^- + 3 -5(l -x)}}. (3.13) 

as expected because in this limit the electron carries all the momentum and the transverse 
width of the impact parameter dependent pdf vanishes . 

For x not equal to 1, TC q (x,b ± ) has nontrivial b ± dependency which comes from the 
(finite) mass terms as well as the A 1 - terms. The q± integration from which we get the scale 
dependency in the above expression (Eqn. I3.13J) cannot be done analytically for this case. 
Here we do not plot the scale dependent TC q (x, b L ). 

Next we consider the helicity flip part of the matrix element : 

<<>) e iP + y~^ P + Aj T | ^(Z|_) 7 +^,(^_) | P, |) = j^x(l - xfi-im^-iA 1 - A 2 ) 



J?r \ 2 2 J (2tt) 

rfV 

((g^) 2 + m 2 (l - x) 2 )((g^ + (1 - x)A^) 2 + m 2 (l - x) 2 ) 
= -^(A'-^). (3.14) 

m is the renormalized mass of the electron. In the limit A 1 - = we get 

OL 

EJx,0) = -x. (3.15) 

7T 

Integration over x gives the Schwinger value for the anomalous magnetic moment of an 
electron in QED : 

/ E q dx = F 2 = — . (3.16) 

JO 2,71 

This result can also be obtained by directly calculating the matrix element of the current 
operator in the light-front framework [l^j ]. 
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Fig. 1: (Color online) (a) £ q (x,b ± ) vs b 1 - for three different values of x; (b) Derivative of 
£ q (x, b- 1 ) with respect to b x as a function of b ± . 

For non-zero A -1 , E q becomes, 

A 2 ) = %f j; *P m _ p) * A±)2 + m2 y (3.17) 

(3 is the Feynman parameter. The scale dependence in this case is suppressed and we can 
integrate over the full | q 1 - | range, from to oo. 
Taking the Fourier transform, we get 



(27T) 2 

AdA 

(2^ 

— AdAE q (x,-A 2 )J (bA). (3.18) 

Z7T JO 



For numerical calculation we use the real integral form of the Bessel function : 

J (bA) = - cos{bAsin6)dO. (3.19) 

7T JO 

If we take a state polarized in the y direction (in the infinite momentum frame) and 
whose center of momentum is at the origin, 

| P + , R L = 0, y) = -^(| P + , R L = 0,l)+i\ P + , R 1 = 0,|)); (3.20) 
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the unpolarized fermion distribution in the impact parameter space gets distorted, 

qy(x, b 1 -) = (P+, R 1 = ± ,y\ O q (x, b 1 ) | P+, R 1 = ± , y) 

= H q (x,b^ + ±-^£ q ( X ,b ± ). (3.21) 

The distortion is directly related to £ q (x, b±). Fig. 1 (a) shows the helicity flip contribution 
S q (x,b ± ) as a function of b for three different values of x. We have plotted for positive 
fe- 1 . £ q (x, b^-) is a smooth function of b 1 - in the range shown and it increases as b 1 - decreases. 
Also, it increases linearly with x, as is clear from Eq. (j3.17j) . We have taken the overall 
normalization ^- = 1 in order to study the qualitative behavior and m — 0.5. S q (x, 
has a maximum at b 1 - = 0. Fig. 1 (b) shows the derivative of S q (x, with respect to b x , 
which gives the distortion of the pdf of an electron at one loop in impact parameter space 
due to transverse polarization. The integrand in this case contains Ji(5A) and it is highly 
oscillatory. However, the integral converges for the b 1 - range shown in the plot. As S q (x, b L ) 
is a smooth function with maximum at b 1 - = its derivative for positive b 1 - is negative. One 
can see that the distortion of the distribution in impact parameter space increases as b L 
decreases and for a given b 1 - the distortion is higher in magnitude for larger values of x. The 
distortion shifts the distribution actually toward negative values of b^. 



IV. GAUGE BOSON DISTRIBUTIONS 



From the definition, 



H g (x,A 2 ) 



87^sP^ 



dy e 



ip+y-x 



(P', 



(4.1) 



we take the state to be a dressed electron as in Eq. ()3.2j) . Contribution in this case comes 
only from the two-particle sector, and we get 



H g (x, A 2 ) = Af 2 J dV^ 2 *(l - x, g ± )^ 2 (l - x, q L + (1 - x)A 

'(i + a-x) 2 



a 

W) 



x 



— h m x 



2„_L 



d 2 q 



x) 2 (A J 



d 2 q J 



d 2 q L | l + (l-o;) 2 
L±L 2 



2x 



d 2 q 1 



(4.2) 



L12 J h\Li2^ 

where = (q ± ) 2 +m 2 x 2 and L 2 = (g _L + (l— x) A ± ) 2 +m 2 x 2 . The scale dependency as before, 
comes from the limits of the q x integrations. Fourier transform of the above expression gives 



J i I i I i L 



0.05 0.1 0.15 0.2 0.25 



(a) 




Fig. 2: (Color online) (a) £ g (x, b ) vs b ± for four different values of x; (b) Derivative of 
£ g (x,b ± ) with respect to b x as a function of b^~. 

the impact parameter dependent gauge boson pdf. Nontrivial b dependence come from the 
mass term and the A x term. The other terms give logarithmic dependence on the scale. In 
the forward limit, 



H g (x,Q) 



a: 
2^ 



1 + (1-X) 2 n Q 2 

/O^— . 

X J jJL z 



(4.3) 



Next we look at the helicity flip part of the matrix element. This can be written as, 



d lL e ^v-^ P + A) T| F+»(^)F+( y -) | P, I) = -^(A 1 - .A 2 ), 

2m 



which gives 



x(l — x) 



(/3(1 - (3)(A^) 2 (1 - x) 2 + m 2 x 2 ) 
If we denote the momentum fraction of the quark as x instead of 1 — x, we get, 

x 2 (l — x) 



In the forward limit, this gives, 



((3(1 - (3)(A^) 2 x 2 + m 2 (l - x) 2 ) 



E g (x, 0) 



a x 

71 1 — X 



(4.4) 



(4.5) 



(4.6) 



(4.7) 
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The second moment of E qt9 (x, 0), J dxxE qj9 (x, 0) gives in units of j- by how much the trans- 
verse center of momentum of the parton q, g is shifted away from the origin in the transversely 
polarized state. When summed over all partons, the transverse center of momentum would 
still be at the origin. Indeed it is easy to check for a dressed electron 

£ dxxE q (x, 0) + J dx(l - x)E g (x, 0) = 0, (4.8) 

which is due to the fact that the anomalous gravitomagnetic moment of the electron has 
to vanish. Note that in the second term, (1 — x) is the momentum fraction of the gauge 
boson. In fact the second moment of E g>g appear in the angular momentum sum rule |2(j] 
and they are related to the orbital angular momentum of the nucleon. In the approach we 
are following, this connection can be seen from the fact that in the light-front gauge, the 
matrix element of the quark orbital angular momentum operator cancels the contribution 
of the gluon helicity and orbital angular momentum for a dressed quark in the helicity sum 



rule 



The impact parameter dependent gauge boson pdf can be defined in the same way as in 
Eq. flUBl): 

= I ^d6E 9 { X ,-tf)e-^ 

= — AdAE g (x,-A 2 )J (bA). (4.9) 
Zix Jo 



Fig. 2 (a) shows S g (x, 6 -1 ) vs b 1 - for three different values of x. E^x^b^) is negative for 
positive b L and has a negative maximum at b 1 - = 0. As before, we took ^- = 1 and m = 0.5. 
Fig. 2 (b) shows the derivative of S g (x, 6 -1 ) with respect to b x as a function of b x . Fig. 3 (a) 
shows £ g (x,b ± ) as a function of x for three different values of b±. Unlike the fermion case, 
S g (x, b- 1 ) for a fixed b L is a complicated function of x and does not increase monotonically as 
x increases. Depending on b- 1 , the maximum of S g (x, b- 1 ) appears for different x (0 < x < 1). 
Again in contrary to the fermionic case, S g (x, b- 1 ) becomes smoother when x is large which 
is expected since for large x (x — 0.7,0.8 in Fig. 2 (b)), gauge boson carries only a small 
(1 — x) fraction of total momentum. Fig. 3 (b) shows the derivative of S g (x, b ) with respect 
to b x as a function of x, which shows the distortion of the pdf in the transverse plane for a 
transversely polarized state. 
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Fig. 3: (Color online) (a) £ g (x, b ) vs x for three different values of b±; (b) Derivative of 
£ g (x, b- 1 ) with respect to b x as a function of x. 

V. SUMMARY 

We investigated the impact parameter dependent parton distributions for a relativistic 
composite system. An ideal framework is based on light-front field theory, where the trans- 
verse boosts behave like Galilean boosts and the longitudinal boost operator produces just 
a scale transformation. We take an effective composite spin 1/2 state, namely an electron 
dressed with a photon in QED. Using the overlap representation of GPDs in terms of light- 
cone wave functions, we obtained the scale dependence of the impact parameter dependent 
pdfs at one loop. We also showed the distortion of the fermion and gauge boson distributions 
in the transverse plane when the state is transversely polarized. Plots show the qualitative 
behaviors of the helicity flip contributions S q (x, b- 1 ) and £ g {x 1 h L ) respectively for the elec- 
tron and the gauge boson and their derivatives with respect to h x which give the distortion 
of the distribution in the transverse plane. For both fermion and gauge boson, for a fixed 
x, the distortions are larger in magnitude for smaller b±. 
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